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Abstract. Bright OH megamaser galaxies, radiating 
1667/1665 MHz lines, could be detected at redshifts from 
z w 1 to 3 in moderate integration times with existing 
radio telescopes. The superluminous FIR galaxies that 
host the megamasers are relatively rare at z w 0, but 
they may have been more common at high redshift, if the 
galaxy merger rate increases steeply with redshift. There- 
fore, blind radio spectroscopic surveys at frequencies of 
400 to 1000 Mhz can form an independent test of the 
galaxy merger rate as a function of time over the redshift 
interval z = 4 to 0.7. 

The redshift range z — 0.17 to ~0.4 will be difficult to 
survey for OH masers, since spectroscopic survey signals 
will be confused with HI emission from normal galaxies at 
redshifts less than 0.3. In fact, the signals from OH masers 
are likely to dominate over 21cm line emission from nor- 
mal galaxies at frequencies below 1200 MHz (i.e. large 
redshifts zhi > 0.18 and zoh > 0.4). Surveyors of nearby 
galaxies in the 21cm line may find that OH masers form 
a contaminant to deep, blind HI surveys for redshift ve- 
locities less than a few hundred kilometers per second. At 
frequencies just above 1420 MHz, sensitive sky surveys 
might detect OH masers, which could be mistaken for a 
population of "infalling, compact High Velocity Clouds" 
but would ultimately be traced to luminous FIR back- 
ground galaxies at z « 0.17 once optical and IR follow-up 
has been performed. 
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1. Introduction 

The detectability of OH megamaser and gigamaser galax- 
ies at cosmological distances has been noted by several 
authors (Baan 1989, Baan et al 1992a, Norman & Braun 
1996, Baan 1997). So far, the sample of known mega- 
masers has been assembled from targeted observations 
of luminous galaxies for which redshifts are known (cf. 
Baan et al 1992b, Martin et al 1989, Stavely-Smith et al 
1992). The receiving systems are tuned to the predicted 
frequencies of the redshifted OH lines whose rest frequen- 
cies fall at 1665 and 1667 MHz. The greatest success 
rate for discovery of megamasers occurs in galaxy sam- 
ples that are selected for having the highest far infrared 
(FIR) 60 \xm luminosity, with the result that ^50% of 
galaxies with L 6 o M > 10 112 /i _2 Lq are strong OH masers 
(h = ^/(lOOkms -1 ). The odds decline to ~5% for sam- 
ples where Lqq^ w 1O 9 L0. There is a striking correlation 
between OH luminosity Loh and Lqq^, with Loh oc Lg 0fl 
(Baan et al 1992a, Baan 1989, Martin et al 1989). 

A picture explaining both the probability of detection 
and the relative strengths of megamaser emission from ob- 
jects of different Lqq^ has the maser strength depending 
on strength of the starburst activity in the galaxy (Baan 
1989, Henkel et al 1991). Major merger and interaction 
events are the strongest stimulants of star formation and 
also produce the most turbulent interstellar gas distribu- 
tions; these conditions provide the strongest IR fluxes for 
pumping the maser levels, as well as the largest covering 
factors of the star forming regions, leading to the greatest 
probability that views from random directions will see OH 
maser activity. Galaxies with milder star formation have 
more placid disks, and maser emission is then observed 
only by observers who view the galaxies nearly edge-on 
to their planes. The strong association between the ultra- 
luminous FIR galaxies and strongly-interacting and merg- 
ing systems has been recently discussed by Clements et al 
(1996). 
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At the present epoch, the superluminous FIR galax- 
ies that host the most luminous masers are relatively rare 
galaxies. This means that radio spectroscopic surveys of 
the sky would need to cover large solid angles before they 
would begin to detect OH masers at random from the lo- 
cal galaxy population. On the other hand, this paper will 
show that, since OH emission can be so very strong, it 
becomes probable that surveys with large spectral band- 
widths and covering large depths are likely to identify cos- 
mologically distant megamasers with moderate integra- 
tion times. 

The most luminous infrared sources known are the ob- 
jects detected at high redshifts, z » 0.4 to 4.7, (Rowan- 
Robinson 1996, and references therein) with Lpm ~ 10 12 
to 10 14 /i _2 Lq. There is general consensus that the merg- 
ing and interaction rate of galaxies was greater in the past 
with merging rate rising in proportion to (1 + z) m . Values 
for m as high as 4 (Carlberg 1992, Lavery et al 1996) are 
mentioned in the literature, although there is some obser- 
vational evidence pointing to m ss 1.2 for z < 1 in the 
HST Medium Deep Survey (Neuschaefer et al 1997) and 
m = 2.8 ± 0.9 from faint galaxy catalogs compiled from 
CFHT surveys (Patton et al 1997). The comoving number 
density of bright QSOs (Mb < —24.5 +51og(/i)) increases 
as (1 + zf from z = to 2 (Hewett et al 1993), and, if the 
QSO phenomenon is tied to galaxy interactions, this steep 
evolution could point to a steep rise in the interaction rate. 

Surveys in the range 400 to 1000 MHz may provide an 
independent measure of the merging rate by determining 
the density of merging galaxies as a function of redshift. 
Since some of the currently most viable models of galaxy 
formation require the bulk of the construction of the larger 
galaxies to occur by merging throughout this range of red- 
shift (z ~ 0.5 to 3), surveys for OH masers will form an 
important test of the principal processes of galaxy forma- 
tion. 

2. Detection rate of cosmologically distant OH 
Masers 

The detection rate for OH megamasers can be estimated 
by combining 1) knowledge of the density of prospective 
FIR luminous host galaxies, 2) the probability that galax- 
ies of each luminosity will be seen as megamaser emitters, 
3) the dependence of megamaser strength on FIR lumi- 
nosity, and 4) the sensitivity of radio telescope receiving 
systems. The steps in the estimation of the detection rate 
are outlined in the following subsections. A parallel line of 
reasoning applies to the detectability of normal galaxies 
in the 21cm line of neutral hydrogen, as discussed in Sect. 
2.5. 

2.1. The FIR luminosity function ^(Lqq^) 

Koranyi & Strauss (1997) have presented a conveniently 
formulated description of the far-infraed luminosity func- 



tion, based on complete samples of galaxies from the IRAS 
60 /im survey. Their paper analyzes the sample in the con- 
text of several cosmological models; here, parameters for 
the "p = 1" case, which is the conventional cosmology, 
are adopted to describe the number density of luminous, 
FIR-selected galaxies. 

The 60 /urn differential luminosity function is given by 

$(L 60 ) =(-?-+ r , f T ) *(£ 6 o) (1) 

with ^(L6o), the cumulative luminosity function, defined 
as 

t(lM) _ c (^)-( 1 + ^)-. 

For the case of conventional cosmology, Koranyi & Strauss 
fit values of a = 0.49, (3 = 1.81 and L* m = 1O 9 68 L . 
The normalization constant, C, is related to their pa- 
rameter m — 0.058 Mpc~ 3 , through the relation m = 
i&[L m i n (z s )], where L m i n (z s ) is the minimum luminosity 
detectable in the sample, whose members are restricted 
to lie at redshifts greater than z s . For Koranyi & Strauss' 
analysis, L min (z s )/L* = 8.8xl0 -3 , which leads to C = 
5.8xl0- 3 Mpc~ 3 (H a = 100 km s" 1 Mpc" 1 ). 

In the Lqq » Lq regime, which will characterize the 
hosts of the brightest OH line emitters, <I>(L6o) may be 
approximated as 

$(i 60 )«C( a + / 3)i* - 1 (i 60 /L* )-^+ 1 ). (3) 

2.2. The effective OH luminosity function Q(Loh) 

The strong correlation between OH and 60/im luminosity, 
Loh oc Lg (cf. Baan et al, 1992a) provides a way to recast 
the 60/im luminosity function for luminous FIR galaxies to 
obtain an OH megamaser luminosity function Q(Loh) as 
shown in Fig. [l]. In this paper, this has been accomplished 
in two complementary ways: one by numerical integration 
of &(L 60 ) to obtain the number of OH sources falling in 
logarithmically spaced bins of Loh, and the second by 
noting that the masers that are detectable at cosmological 
distances inhabit host galaxies that are far more luminous 
than Lg and then using the approximation in Eq.(^|) to 
derive an analytic expression. 

The advantages of computing Q(Loh) numerically are 
that it (1) simplifies inclusion of the spread of OH lumi- 
nosities in a band of width ^10 ±a7 about the Loh oc L\ q 
relation (Baan et al 1992a), and (2) allows for varying the 
probability f m that luminous FIR galaxies will appear to 
be megamasers as a function of Lqq. Here, the range of 
FIR luminosities that are considered to contribute OH 
megamasers is conservatively restricted to the range where 
strong OH emission has actually been detected (cf. Fig. 2 
of Baan et al 1992a); adopting H a = 100 km s _1 Mpc _1 , 
this range is 5xl0 9 < L 6Q < 2xl0 12 i Q . Within this 
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Fig. 1. Estimate for the effective OH megamaser luminosity 
function at z ~ 0. The vertical axis gives the apparent number 
of observed per Gpc 3 between 10" 05 L O h and 10+ 05 Loh, 
with log(Loff) in solar luminosities plotted along the horizon- 
tal axis. The approximate FIR luminosity, Lao of these galax- 
ies is labeled along the top border. With a observational de- 
tection limit of 1 mjy in the OH lines, only a portion of the 
high luminosity side of the function would be detectable at 
high redshifts, as indicated by the labeled cut-off luminosi- 
ties. The portion detectable at z = 1 is drawn as a bold 
line. (H = 100 km s -1 Mpc^ 1 ). The dashed line is the ap- 
proximation for O(Lqh) from Eq.(5) with f m = 1.0. 



range, f m takes on non-zero values as follows, f m — 0.05 



for 5xl0 9 < L 60 < 2xlO iu L Q , f m = 0.2 for 2xlO io < 
L 6Q < 2xlO n L , and f m = 0.5 for 2xlO n < L 60 < 
2xl0 12 L Q (from Baan 1989, Baan et al 1992b). Since 
these probabilities are likely to be a result of viewing an- 
gle, the derived luminosity function should be considered 
as an "effective OH luminosity function" that provides 
the probability of detection for randomly oriented sources, 
rather than an an accurate accounting of the number of 
OH emitting galaxies. The constant of proportionality C\ 



is 4.5x10" 21 Lq 1 
can be defined to be L* OH = C\L*^ = 0.10X 



in Lqh — Ciig 



and a fiducial L 



OH 



Figure [l] shows the numerically computed luminosity 
function for OH megamasers in logarithmically spaced 
bins of 0.1 Dex. For reference, the 60 /im luminosity 
(^60 = V Loh /Ci) is indicated on the top axis. Since 
the bright OH megamasers that are likely to be useful in 
cosmological studies lie well above L* OH , the approxima- 
tion in Eq. (j|) can be used to verify the numerical estimate 



in the range of greatest interest: 

0{L O H)dL O H = f m $[L m {L H)] 

{a + (3) 



dL 



(ii! 



dL 



OH 



dL 



OH 



(4) 



OH 



J OH 



dLpH 
1 



OH 



Since it is convenient to plot the number of objects per 
logarithmically spaced bins of luminosity, the number of 
objects falling in bins of ALqh = 0.1 Dex centered on 
Lqh can be obtained from 
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This function, with f m set equal to 1, is plotted in Fig. [j] 
for comparison with the result of the numerical integra- 
tion. 

In a related analysis, Baan (1997) presents the 60/i 
FIR luminosity function for OH megamaser sources as a 
function of Lqq^ for direct comparison with the luminosity 
function for broader population of all FIR galaxies. Baan 
(1997) further bins the fraction f m in 0.5 Dex rather than 
the 1 Dex divisions used in the approximation here. The 
finer binning would lead to less abrupt inflections in the 
OH luminosity function derived here (cf. Fig. |l|). Overall, 
there is good agreement between the spatial densities of 
galaxies between Baan (1997) and the present work, al- 
though the fraction of OH galaxies showing megamaser 
emission presented by Baan (1997) appears to be roughly 
a factor of two lower in each luminosity range than in the 
earlier publications. This factor of two uncertainty can be 
propagated through the subsequent estimates for volume 
densities and densities per solid angle. 

2.3. Profile width and detectability at large redshift 

The integral line flux observed from a cosmologically dis- 
tant emitter at redshift z is (Wieringa et al 1992) 



Loh 
L 



(6) 



OH 



Ho 
4c 



[n oZ + (n„ - 2)yi + n oZ - 1)] 2 



where is the luminosity distance, H D is the Hub- 
ble constant and fl is the cosmological mass density 
normalized to the critical density. If the line of width 
AV km s _1 is emitted at rest frequency v and rcdshiftcd 
to v — v /{\ + z), the observed frequency width of the 
line will be Av w v AV/ c, providing an average flux den- 
sity over the line profile of S u w S/ Au cx S(l + z). Thus, 
the line width is squeezed in frequency as a line of velocity 
width AV^ is observed toward higher redshifts, leading to a 
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partial compensation for the reduction in the integral line 
flux by the inverse square law (Eq.(6)), The range of line 
widths observed for OH megamaser galaxies is illustrated 
in Fig. ||, where the line width is plotted as a function of 
maser luminosity. The vast majority of the line profiles fall 
in the range 30 to 300 km s _1 . Only the most luminous 
source known shows a much broader profile (Baan et al, 
1992a). 





Fig. 2. Velocity widths of OH megamaser profiles as a function 
of luminosity in the 1667 OH line. Data from Staveley-Smith 
et al 1992, Martin et al 1989, Baan et al 1992a, Baan et al 
1992b. 



The sensitivity of a radio telescope is generally speci- 
fied (cf. Crane & Napier 1989) in terms of its receiver noise 
level as parameterized by a system temperature, T sys , and 
the antenna sensitivity, K a = lO~ 26 e ap A/2k, which is a 



Fig. 3. Minimum Loh detectable at redshift z in an observa- 
tion with detection threshold Sdet ■ Families of curves are drawn 
for Sdet = 0.2, 1, 5, and 20 mjy for cosmological models with 
tt = 0.01, 0.2 and 1. Models assume H a = 100 km s" 1 Mpc" 1 
and OH profiles having width 100 km s - . 



with adequate spectral resolution for the identification of 
100 km s _1 width signals. The minimum detectable OH 
line luminosity as a function of redshift is plotted in Fig. |^ 
for a range of cosmological models and observational sen- 
sitivity. 



relation between flux density, S v in Jy and antenna tem- 
perature T a = K a S u K. The quantities e ap , A, and k are 2 -4- Density of megamasers per solid angle 



c ap> 

the aperture efficiency, the antenna collecting area and 
the Boltzmann constant in mks units, respectively. The 
noise level attained by observing a single polarization with 
a spectral resolution of Av and an integration time r is 
ajy — T sys K c /(K a y/ Av r) where K c is a factor with val- 
ues typically between 1 and 1.4 depending on the details 
of the digital correlation spectrometer used to obtain the 
spectra. 

As an example of sensitivities that can be attained at 
present or in the near future, consider the specifications 
and noise levels for typical observations with the Wester- 
bork Synthesis Radio Telescope. For this system, receivers 
at 600 and 1400 MHz will have T sys of 45 and 25 K re- 
spectively. The effective K a of combining fourteen 25 m 
telescopes is K a w 1.4 K Jy -1 . For velocity widths of 100 
km s^ 1 , these receiving systems will reach the 1 mjy (5er) 
level in roughly one to three 12 hour long integrations. 
The WSRT telescope will soon have the further benefit of 
a spectrometer capable of simultaneous observation in the 
synthesis mapping mode of an 80 MHz bandwidth at fre- 
quencies below 1200 MHz and 160 MHz above 1200 MHz 



Estimates for the number of detectable megamasers in a 
blind survey of a large area of sky can be straightforwardly 
obtained by computing the comoving volume at redshift 
z contained in the solid angle dQ and depth dz (Wieringa 
et al 1992), 



dV = 4 



c 

W 



[n z + (n -2)(yi + n z-i) 
ni(i + z) 3 ^i + n z 



-dz dfl. 



and then filling the volume with OH emitters according to 
the density prescribed by Fig. 0. A count of the number of 
objects whose observed flux exceeds a detection threshold 
can be made and plotted as shown for three cosmological 
models in Fig. [|. The detection level was chosen to be 1 
mjy for 100 km s _1 , consistent with feasible integration 
times with existing radio telescopes. In Fig. Q the horizon- 
tal axis is chosen to be observed frequency and the counts 
are binned according to expected detection rate per band- 
width of radio spectrum. These arc natural units for use in 
planning surveys. In computing this diagram, the assump- 
tion was made that there is no evolution in the comoving 
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number densities or luminosities of the OH megamasers 
with rcdshift. Curves for showing the detection rate for 
the confusing signals from normal galaxies emitting the 
21cm hydrogen line are shown for comparison. Computa- 
tions of the detection rate for a range of sensitivities is 
given in Fig. [5]. 
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Fig. 4. Detection rate of signals from OH megamasers and 
HI in normal galaxies as a function of frequency. Computed 
for line widths of 100 km s _1 and survey sensitivity of 1 mjy. 
The OH redshift for the OH lines is indicated at the top of the 
figure. The relatively narrow peak in detection rate for neutral 
hydrogen is labeled "HI," while "OH" refers to the broader, 
low-level function of redshift. Three cosmological models are 
shown: n o = 0.01, 0.2 and 1. 



2.5. Confusion with 21cm line emission from normal 
galaxies 

The expected number of detections of normal galaxies in 
the 21cm line of neutral hydrogen can be computed in a 
similar manner to that for OH megamasers. The HI lu- 
minosity function 4>(Mhi) observed for nearby galaxies 
(Zwaan et al 1997) can be represented as a function of HI 
mass, Mhi 



4>(M HI ) dM m = 0* (M HI /M* Hi y 



-M HI /M' HI dMsi ( 7 ) 



with <p* = 0.014 Mpc" 3 , 7 = 1.2 and M* HI = 10 9 55 M Q . 
The luminosity in the 21cm line from an optically thin 
cloud of neutral hydrogen with HI mass M^i is L^i/ Lq — 
6.2xlO- 9 M HI /M Q . 

Curves for detection rates in the hydrogen line are plot- 
ted in Figs. [| and || for comparison with the OH line 



strengths. The exponential cutoff to the HI luminosity 
function causes the HI detection rate to cutoff very hard 
around z — 0.2 for 1 mjy sensitivity. The power law func- 
tion for the OH masers permits the bright end of the 
OH maser population to dominate the detection rate at 
frequencies below ^1200 MHz for reasonable integration 
times with existing telescopes. 

Not only are the typical velocity widths of the OH pro- 
files comparable to velocity spreads measured in the 21 cm 
line for low mass and face-on galaxies, but also the pair 
of strong OH lines at 1665.4 and 1667.4 MHz could pro- 
duce a spectrum that, at low signal to noise ratio, could 
imitate a double horned galaxy profile with a splitting of 
^360 km s -1 , although the horns are likely to be asym- 
metric due to the usual weakness of the 1665 line relative 
to the 1667 line. 

Sifting the OH masers from the HI lines from nor- 
mal galaxies of low redshift would be straightforward. 
High spatial resolution radio observations would show OH 
masers to have high brightness temperature and would 
identify the sources with compact optical sources for which 
emission lines would provide an unambiguous redshift. 
OH megamasers are strong emitters of 60 /im infrared, 
but sources at z 0.17 (or greater), selected for OH line 
strength around 1 mjy, will have 60 /xm flux densities ap- 
proximately equal or less than the 0.2 Jy sensitivity level 
of the IRAS 60 /im Faint Source Survey (Moshir et al 
1992). Since the infrared spectrum of the OH megamasers 
falls steeply from 60 [iva. to 25 /im, suitable megamaser 
host galaxies at moderate to high redshift may become in- 
creasingly difficult to detect at infrared wavelengths as the 
60 fj,m flux is redshifted toward still longer wavelengths. 

3. Redshift dependent merging rate 

The importance of merging small galaxies to build the 
large galaxies that we observe at the present is a topic 
of considerable controversy, since the question addresses 
the foundations for theories of galaxy formation. In this 
section, it is demonstrated that surveys sensitive to OH 
masers at z w 2 may be able to measure the number den- 
sity of interacting and merging systems at high redshift 
at a time when the merging rate is thought to peak. Ob- 
servations using radio telescopes with sufficient sensitiv- 
ity and sufficient spectral bandwidth can not only deter- 
mine the statistics of the evolving density of OH masers 
with time, but they will also select systems of interest at 
high rcdshift, which can be further studied with a variety 
of techniques to determine physical conditions, dynamics 
and effects of local environment. 

For the purposes of estimating the effect of an in- 
creased merging rate at early times, the model explored in 
Fig. fl| and the top panel of Fig. || can be modified by in- 
creasing the number density of OH megamasers in propor- 
tion to (l + z) m , in keeping with the idea that, locally, OH 
megamasers are associated with luminous FIR galaxies, 
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Fig. 5. Detection rate as a function of frequency for two dif- 
ferent redshift dependences of the merging rate. Top Panel 
Constant merging rate. Bottom Panel Merging rate propor- 
tional to (1 + z) 4,5 for 2 < 2.5 but constant for z > 2.5. Detec- 
tion rate is computed in a cosmological model with fl = 0.2 
for four sensitivity levels (0.2, 1, 5, and 20 mjy) for signals of 
velocity width 100 kms -1 . The lowest sensitivity survey (20 
mjy) corresponds to the most heavily shaded area, and the 
highest sensitivity (0.2 mjy) is unshaded. 



which in turn are found to occur in systems with merging 
and interacting morphologies (Clements et al 1996). For 
comparison, an extreme example, m = 4.5 (Lavery et al 
1996) is shown in Fig. ||. The significance of this detection 
rate is that an existing telescope, such as the Westerbork 
Synthesis Radio Telescope, which will reach a noise level 
below 1 mjy level in a 12 hour synthesis observation and 
can survey a field of one square degree per exposure at 600 
MHz (z h ps 2), would be likely to detect a few objects per 
field, when it is equipped with a spectrometer capable of 
80 MHz (or more) bandwidth. 

The luminous quasar population shows a still stronger 
evolution than is typically suggested for normal galaxies. 
Since the most luminous OH masers are associated with 
the most luminous galaxies in the nearby universe, it is 
reasonable to speculate about a possible relation between 
hosts for OH maser and quasars and the possibility that 



both may be driven by mergers and interactions. Figure |^ 
shows recent observational constraints on the comoving 
number density of bright QSOs (M B < -24.5 +5 log(/i)), 
as presented by Hewett et al (1993) for z < 2 and Schmidt 
et al (1995) for z > 2.5; these studies provide evidence 
that the comoving density peaked at z ss 2 to 3. Rep- 
resentative curves for describing the increase in density 
</>_24.5 with redshift are drawn in the figure and indicate 
that 0_24.5 is roughly proportional to (1 + z) 6 to redshift 
as high as z sa 2. The QSO evolution is often described 
as an evolution in luminosity; Boyle et al (1988) quan- 
tify the QSO population as having a luminosity function 
cf> qso (L/L*) cx </> (L/L*)- 3 - 7 with L* = L D (1 + z) 3 - 2 . In- 
tegrating to obtain the number of objects per comoving 
volume brighter than L\ gives 

"^'■jC>(e)'(e) (8) 

which has a proportionality that is consistent with the 
evolution for z < 2 shown in Fig. ^| and in Schmidt et al 
(1995). 

The intention behind Fig. || is to show that at least 
one population of luminous object does evolve very steeply 
with redshift. An association between AGN activity and 
galaxy merging or interaction remains to be definitively 
established, although there is evidence that "warm," ul- 
traluminous infrared galaxies may be an evolutionary step 
in the evolution of optically selected QSOs (cf. Surace et 
al 1998). Interestingly, the z « number density of ultra- 
luminous FIR galaxies is comparable to the z w density 
of luminous QSOs, as illustrated in the figure. Further- 
more, Saunders et al (1990) deduced a density evolution 
oc(l + z) 6 7 for their 60/zm sample, although subsequent 
study by Ashby et al (1996) did not find the tail to high 
redshift expected if this were indeed the case. The evolu- 
tion rate oc(l + z) 4 ' 5 explored in the lower panel of Fig. || 
is fairly mild in comparison to that of the QSOs. 

4. Contamination of blind HI 21cm line surveys 

Figures. |] and ^| demonstrate that blind spectroscopic sur- 
veys are more likely to detect distant OH megamasers than 
HI from normal galaxies in the frequency bands below 
~1200MHz, assuming that both megamasers and normal 
galaxies have similar properties to those that are observed 
nearby. Above the rest frequency of the hydrogen line at 
vjij = 1420.4 MHz and in a small range of frequencies just 
below Vju there is a possibility for confusing OH emission 
at z OH ~ 1667/1420 - 1 = 0.17 with weak HI signals 
from nearby dwarf galaxies and HI clouds. To explore this 
possibility in more detail, Fig. shows a a plot similar to 
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Fig. 6. Comoving density of bright QSOs (Mb < —24.5 
+51og(/i)) as a function of redshift. Points at z < 2.5 are 
derived from Hewett et al (1993) (converted to H = 100 
km s _1 Mps -1 ); points at z > 2.5 are adapted from Schmidt 
et al (1995). Three curves are drawn to indicate the evolu- 
tionary trend through the low redshift regime with the depen- 
dence oc (1 + z) n with n = 2, 4 and 6.5. A horizontal dotted 
line is drawn to indicate a constant comoving density equal 
to the current density (see Fig. [j]) of luminous OH sources 
(Loh > lO 41 fr~ 2 L and L 6 o M > 1O 12 ' 2 /i~ 2 L ) that are bright 
enough to be detected at z = 3 with 1 mjy sensitivity level. 



Fig. |5| with an expanded scale around 1420 MHz. The ver- 
tical axes have been adjusted to more convenient units to 
address this question. The possibility for confusion might 
arise in "blind surveys" for dwarf galaxies in the HI line, 
where signals at noise levels of a few a are selected as can- 
didate dwarfs. At these sensitivity levels, only a narrow 
1667 MHz might be detected, mimicing a narrow HI pro- 
file of a few 10's of km s _1 velocity width. Clearly more 
sensitive followup in the radio line, as well as far infrared 
and optical observations, would quickly show the distinc- 
tion. 

Some questions remain about the slope of the luminos- 
ity function for galaxies at the faint end, and this uncer- 
tainty propagates into the faint galaxy tail of the Hi-mass 
function as well. Recent HI surveys have produced val- 
ues for 7 (as defined in Eq.(^|)) in the range 1.2 to 1.7 
(Zwaan et al 1997, Schneider 1997), although there is also 
observational evidence for a sharp decline in the number 
density of HI rich objects at faint HI luminosity (Hoffman 
et al 1992), at least in the Virgo cluster. Figure implies 
that surveys with sensitivities of ~1 mjy or better may be 
likely to detect OH masers (at zqh ~ 0.17) in comparable 



EC 

o 
m 

CD 

a 

CO 

TS 

CO 

CD 



CD 
ft 



CD 

6 

z 

an 
o 




1415 



1420 



1425 



2000 



1000 













y = 1.3 ~ 



















































-1 

-2 



N 

EC 
ES 



CD 

a 



cm 

<D 

a 

cr 

w 

u 

CD 

ft 
u 

CD 



3 

55 



60 
O 



1410 1415 1420 

Frequency [MHz] 



1425 



Fig. 7. Detection rate of HI and OH signals near the rest 
frequency of neutral hydrogen. Two cases are illustrated for 
the faint end slope of the HI mass function, 7 = 1 and 1.3. 
Detection thresholds are set at 0.2, 1, 5 and 20 mjy. 



numbers to the rate at which they detect weak HI signals 
in the velocity range vri < 500 km s _1 . 

This regime can be explored in more detail by esti- 
mating the form of the HI mass function that the OH 
masers would imitate. Figure |^ is an aid to visualize the 
problem. It shows a volume dVoH at zoh that is mapped 
into the volume dVjji at z due to the overlap that occurs 
in observed frequency for the two populations. For sim- 
plicity in this estimate, approximations for the luminos- 
ity distance di ~ zohc/H q and for the comoving volume 
dVoH ~ ZQ H dzondfl will be assumed, instead of the more 
rigorous forms used earlier in this paper in the discussion 
of OH masers at cosmological distances. These produce er- 
rors of less than 10 percent in di and less than 40 percent 
in comoving volume. 

The incremental frequency band dv observed at the 
frequency v is related to increments in redshift for the 
OH and HI lines: 



dv = dzoH = dzfji 



(9) 



implying that dzoH — (yon /vui)dzHi- The flux S mea- 
sured in the line emission from the OH megamaser of lu- 
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Fig. 8. The volume dVoH « ZQ H dzondQ containing OH 
masers that is mapped into volume dVui = z 2 dz H id£l in HI 
surveys. 

minosity Loh is falsely interpreted as HI emission from 
an object of Lhi- 



S = 



Lqh 



-iHI 



Mz 0H c/H )2 Att(zc/H Y 



(10) 



so that Lhi = (z/ ' zoh) 2 Loh- Since the redshift range 
being explored for imitation HI signals is small, z < 0.01, 
the OH masers are drawn from a thin shell, allowing zoh 
to be considered constant at zoh ~ 0.17 in what follows. 

The goal is to derive an approximation for the imita- 
tion HI mass function <j> im provided by distant OH masers. 
This can be accomplished by evaluating the number of sig- 
nals that appear to come from the volume dVui from emit- 
ters of strength L H i but which were actually produced 
in dVoH from emitters of L Q h- The relation between 
the luminosity functions then is Q(Lon)dLoHdVoH = 
<t>im{LHi)dL,HidVHi- Solving for (j> im , substituting from 
the above expressions as necessary and then including the 
approximate relation for the OH luminosity function from 
Eq.(4), which is valid here since all the OH signals de- 
tected will originate in emitters much brighter than L* OH : 



- m 

= e((^) 2 ^)(^ 



(11) 



VQH 
VHI 



C(a + /3)f m v OH ( L HI \ (ZO H\ 2 - a -P 
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This result can then be recast as a dependence on Mhi- 



tiMfidMm = Cfr 



a + P ( L* OH 



2 V L H i 
(MhiY 



VQH 



vhi \M* m ) 
I z \ dMni 



\ZOH 



M 



(12) 



After substitution for the constants in Eq.(12), 
the functional dependence 



t(M HI )dM H i oc 



M 



HI 



M* 
m HI 



-2.15 



ZOH 



0.3 



dM, 



HI 



M 



has 



(13) 



HI 



producing a form with that mimics a diverging low lumi- 
nosity tail, since (a + [3 + 2)/2 = 2.15 is greater than 2. 
The dependence on z does not affect the power law shape 
of the imitation mass function but does weakly affect the 
normalization that would be inferred. 

The number of imitation HI masses per decade of Mhi 
becomes 



2.3x10" 



M 



HI 



M 



HI 



-1.15 



zoh 



0.3 



Mpc~ 



(14) 



HI 



This estimate implies that a survey sensitive to 1 mJy 
signals, corresponding to Mhi between 10 5 and 10 6 M Q at 
velocities of 300 to 500 km s" 1 , would be unlikely to detect 
more than one imitator per cubic megaparsec. Therefore, 
with the current set of parameters describing the OH lu- 
minosity function, significant contamination of HI surveys 
at low redshift would only occur if the HI mass function 
turns out to be as flat as 7 w 1 or drops off at low Mhi- 
Even substantial evolution in the merging rate from z = 
to 0.17 is unlikely to raise this density estimate by more 
than a factor of two, but some caution is in order until our 
understanding of the properties of OH masers improves. 
The likelihood for confusion might rise significantly if sur- 
veys were conducted at a 0.2 mJy sensitivity. 

5. Conclusion 

The density of OH megamaser galaxies in the sky may 
be high enough that radio spectroscopic surveys will be 
effective tools at identifying them. The time evolution of 
the megamaser sources could then be monitored by a series 
of surveys performed at a range of radio frequencies (i.e. 
different rcdshifts). 

OH megamasers may form a source of confusion in sur- 
veys designed to detect neutral hydrogen in normal galax- 
ies through their 21 cm emission. High resolution radio 
mapping, optical spectroscopy and far infrared detection 
would remove the ambiguity. 
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